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Abstract: Analysis of environmental significance and hydrochemical characteristics of river water in 
mountainous regions is vital for ensuring water security. In this study, we collected a total of 164 water 
samples in the western region of the Altay Mountains, China, in 2021. We used principal component 
analysis and enrichment factor analysis to examine the chemical properties and spatiotemporal variations 
of major ions (including F-, CF, NOs, S047, Lit, Nat, NH4*, K*, Mg?*, and Ca?*) present in river water, 
as well as to identify the factors influencing these variations. Additionally, we assessed the suitability of 
river water for drinking and irrigation purposes based on the total dissolved solids, soluble sodium 
percentage, sodium adsorption ratio, and total hardness. Results revealed that river water had an alkaline 
aquatic environment with a mean pH value of 8.00. The mean ion concentration was ranked as follows: 
Ca2*>SOu?->Nat>NO3->Me?*>Kt>Cl>F->NHy4t>Li*. Ca?*, SO4?, Na*, and NOs" occupied 83% of 
the total ion concentration. In addition, compared with other seasons, the spatial variation of the ion 
concentration in spring was obvious. An analysis of the sources of major ions revealed that these ions 
originated mainly from carbonate dissolution and silicate weathering. The recharge impact of precipitation 
and snowmelt merely influenced the concentration of Cl, NOs, SO4?, Ca?*, and Nat. Overall, river 


water was in pristine condition in terms of quality and was suitable for both irrigation and drinking. This 
study provides a scientific basis for sustainable management of water quality in rivers of the Altay 
Mountains. 
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1 Introduction 


Rivers, which support vital ecosystems, serve as crucial surface water resources that supply 
drinking water and irrigation water to the public (United Nations, 2019). In recent decades, the 
quality of river water has faced mounting pressures due to the impacts of human activities and 
climate-related stressors, leading to noticeable water environmental issues worldwide (Wen et al., 
2015; He et al., 2020; Yapiyev et al., 2021). Consequently, there is an urgent requirement for 
continuous river monitoring to ensure sustainable management of water quality. The chemical 
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composition of river water is an essential indicator of water environmental quality, regional 
environmental characteristics, and the distribution and transformation of water components (Ye et 
al., 2010; Mapoma et al., 2016; Niu et al., 2017). Monitoring the chemical conditions of river 
water provides valuable insights into regional chemical weathering and key processes that govern 
hydrochemistry (Wang et al., 2015; Iqbal et al., 2018). Additionally, analyzing the chemical 
properties of river water aids in identifying the sources of major chemical constituents and 
understanding the evolution of geochemical solutes (Li et al., 2017; Pant et al., 2017). 

Studies on the chemical status of runoff water have revealed the major contributors affecting 
river water quality, including atmospheric inputs, rock weathering, and human activities (Li et al., 
2019). The solutes in river basins (particularly in the glacial catchments) predominantly originate 
from crustal mineral weathering (Li et al., 2022). Several qualitative and quantitative methods are 
used for assessing the physicochemical parameters and sources of solutes in river water. For 
instance, qualitative methods, such as Piper histograms (Piper, 1944), Gibbs scatter plots (Gibbs, 
1970), factor analysis (Vermette et al., 1988; Okay et al., 2002; Li et al., 2014), and ionic ratios 
(Villegas et al., 2013), have been applied to identify the runoff water chemistry and mineral 
weathering. Quantitative methods, on the other hand, have been used to quantify the contribution 
of different solute sources based on chemical concentrations. For example, the forward model 
(Galy and France-Lanord, 1999) was utilized to determine the contribution of riverine solutes and 
could accurately assess the differences in the status of crustal mineral weathering (Yu et al., 
2021). 

Numerous previous studies have investigated the hydrochemistry of prominent rivers 
worldwide, including the Nile River (Dekov et al., 1997), Amazon River (Stallard and Edmond, 
1983), and Ganges River—Yarlung Zangbo River (Galy and France-Lanord, 1999). In China, 
researches predominantly focused on hydrochemical characteristics of rivers originating from the 
Tibet Plateau (Li et al., 2019), Pamir Plateau (Wu et al., 2020), Qilian Mountains (Li et al., 2014), 
Tianshan Mountains (Yapiyev et al., 2021), and Altay Mountains (Liu et al., 2021). These 
researches have provided crucial insights into regional environmental changes, such as variations 
in the spatial and temporal distribution of river water chemistry (Jiang et al., 2015), as well as the 
sources of solutes in river water and the factors governing them (Li et al., 2019; Dong et al., 
2022). Evaluating the hydrochemistry of rivers in these regions is vital for effective and 
sustainable management of water quality. 

The western region of the Altay Mountains represents the northernmost concentration of 
modern glaciers in China (Wang et al., 2015). It not only serves as a representative arid region 
influenced by westerly circulation but also functions as a significant water source conservation 
area in northern Xinjiang Uygur Autonomous Region of China. In recent years, researchers 
have increasingly focused on water resources and environmental issues in the Altay Mountains, 
including river water chemistry (Liu et al., 2021), river water quantity (Zhang et al., 2010; Li et 
al., 2018), and soil-related concerns (Goenster-Jordan et al., 2021). However, studies on river 
water chemistry have typically relied on single-time-period sampling results to interpret 
concentration patterns, neglecting seasonal and annual variations (Liu et al., 2021). 
Furthermore, there is a lack of information on spatial patterns, such as the elevation distribution 
of ion concentration, due to scattered data and sampling discontinuities. The seasonal evolution 
and controlling factors of major ions in river water have not been systematically evaluated. To 
explore the hydrochemical characteristics and controlling factors in this region, researchers 
established the Altay Observation and Research Station for Cryospheric Science and 
Sustainable Development in the Altay Mountains in 2016. The primary objectives of this study 
were to: (1) analyze the spatiotemporal variability of major ions in river water, (2) determine 
the chemical weathering processes and the mechanisms underlying the seasonal evolution of 
river water in the study area, and (3) assess the suitability of river water for drinking and 
irrigation purposes. 
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2 Materials and methods 


2.1 Study area 


The western region of the Altay Mountains is mainly the origin area of the Burgin River 
(86°48'-88°36'E, 47°42'-49°12’N) and Haba River (85°31'-87°09'E, 47°38'—49°09'N; Fig. 1). 
The Burgin River is the largest tributary of the Ertix River, and the Burqin River Basin has a 
drainage area of approximately 8422.0 km’, with an elevation of 470-4300 m. The glaciers in the 
Youyi Peak region (48°40'-49°10’N, 87°36'—87°53'E) supply water to the Burqin River Basin. 
The Haba River is the second largest tributary of the Ertix River. The drainage area of the Haba 
River Basin is approximately 7224.0 km?. In 2020, the Altay Mountains contained 1927 glaciers 
with a total area of 1096.0 km? (Chang et al., 2022). The glaciers with an area of 0.1—-0.5 km? 
(815 glaciers) account for 42% of the total number of glaciers (Chang et al., 2022). The Altay 
Mountains above 3000 m is mostly covered by glaciers and seasonal snow, which begin melting 
in June when the air temperatures are above 0.0°C (Qin et al., 2020). The landscapes are 
characterized by glaciers, snow, permafrost, forestlands, and grasslands from high to low 
elevations. 


DEM (m) 
mam High: 8806 
- Low: 443 
—— Air masses 
>» Sampling site 
—— River 
— i HE Glacier 


Fig. 1 Overview of the study area and distribution of the sampling sites. The Digital Elevation Model (DEM) is 
derived from www.earthdata.nasa. gov. 


The entire region has a typical temperate continental climate (Liu et al., 2021). The local 
climate may also be influenced by westerlies and polar air masses. Besides, the annual average air 
temperature is 5.0°C, with a maximum temperature of 39.8°C occurring in July and a minimum 
temperature of —27.9°C occurring in December. The average annual precipitation is 159 mm. 
Precipitation occurs mainly due to water vapor transported by the westerly airflow throughout the 
year (Wang et al., 2015). These climatic characteristics are recorded by the meteorological station 
(86°87'E, 47°71'N; 483 m), which is situated about 5 km northeast of the confluence of the 
Burqin River and Ertix River. In addition to precipitation, glacier meltwater also supplies the 
region. The meltwater accounts for 45%-50% of all water sources (Lanzhou Institute of 
Glaciology and Geocryology, 1982). 
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2.2 Sample collection and measurement 


In order to investigate hydrochemical characteristics of rivers in the western region of the Altay 
Mountains, we collected 164 samples from this region during January—December in 2021. Field 
sampling was conducted every 2 h from the 20" to the 21% of each month. The sampling sites are 
depicted in Figure 1. Firstly, empty sample bottles were washed by using a dilute HCI solution. 
After sampling, the collected samples were filtered using a 0.45-um Millipore filter membrane. 
Subsequently, each sample was poured into two pre-washed polyethylene bottles. One of the 
bottles, which was used for cation analysis, was acidified using HNO; to a pH value below 2.00. 
The other bottle was used for anion and isotopic analyses. All samples were stored at a 
temperature of approximately 4.0°C before and during transportation to the laboratory for 
analyses in December 2021. Each sample was naturally thawed at normal room temperature 
(about 15.0°C) before being determined in the laboratory. 

The hydrochemical analysis was conducted at the State Key Laboratory of Cryospheric Science 
of the Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences. The 
pH was measured by a pH meter (PHSJ-4A, INESA Scientific Instrument Co., Ltd., Shanghai, 
China). We determined the total dissolved solids (TDS) and electrical conductivity (EC) based on 
a conductivity meter (DDSJ-308A, INESA Scientific Instrument Co., Ltd., Shanghai, China). The 
concentrations of cations (including Lit, Na*, NH4*, K*, Mg?*, and Ca*) were measured using an 
atomic absorption spectrometer (PE2380, Perkin Elmer, Massachusetts, USA), while those of 
anions (including F`, Cl, NO7, and SO4%) were determined using an ion chromatography 
apparatus (Dionex100, Thermo Fisher Scientific, Massachusetts, USA). To ensure that the error 
between the analysis results of the water samples and those of the standard sample was less than 
5%, we calibrated the apparatus after every 20 sample measurements. 


2.3 Statistical methods 


The principal component analysis (PCA) is a robust multivariate analysis method employed to 
identify the significant parameter variations within the data (Shrestha and Kazama, 2007; Wu et 
al., 2020). In this study, the PCA was utilized to condense several interconnected major ions 
(i.e, F, Cl, NO7, S047, Lit, Na*, NH4*, K*, Mg?*, and Ca?*) into fewer uncorrelated 
variables, thereby facilitating subsequent analyses. The PCA algorithm was carried out in SPSS 
19.0 software. 

To assess the role of mineral weathering as the source of ions in this region, we calculated the 
enrichment factor (EF) using Ca? and Na” as the reference material. Ca?* and Na* serve as the 
reference materials for soil and sea, respectively. The EF was calculated as follows: 


EF iI =[ X / Ca” fivater sample I[X / Ca% Joi? (1) 


EF ea = [x / Wa ocru /Lx / Na* hea ý (2) 


where EFgoi is the enrichment factor for soil; X is the concentration of a specific ion (mg/L); 
[X/Ca**]water sample is the ratio of the concentration of a specific ion to the concentration of Ca?* in 
the water samples; [X/Ca?*]soi is the ratio of the concentration of a specific ion to the 
concentration of Ca** in soil (Okay et al., 2002); EFsea is the enrichment factor for sea; 
[X/Na™ water sample is the ratio of the concentration of a specific ion to the concentration of Na* in 
the water samples; [X/Na*]sea is the ratio of the concentration of a specific ion to the 
concentration of Na” in sea (Keene et al., 1986); and Na” (mg/L) and Ca* (mg/L) represent the 
concentration of Na* and Ca?*, respectively. 


2.4 Suitability indices for drinking water and irrigation water 


High concentration of Na” in irrigation water can lead to sodium hazards. Na* replaces both Ca?* 
and Mg”, resulting in lower permeability and solid hardening (Shaki and Adeloye, 2006). In this 
study, we assessed the suitability of the water samples for irrigation based on the soluble sodium 
percentage (SSP; %) and sodium adsorption ratio (SAR) (Wu et al., 2020), which were calculated 
using the following equations: 
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; 

SSP = bi - —x100%, (3) 
Na* +K* +Ca** +Mg** 

SAR = Na*/,[(Ca?* +Mg”*)/2. (4) 


where K* (mg/L) and Mg” (mg/L) represent the concentration of K* and Mg”, respectively. 

In addition, the total hardness (TH; mg/L), which refers to the combination of carbonate and 
non-carbonate hardness (Gao et al., 2017), was determined using the following equation: 

TH =| Mg” /24+ Ca? /40 |x100. (5) 

On the basis of the TH values, we classified the water samples into the following five classes: 
very soft water (TH<75 mg/L), soft water (75 mg/L<TH<150 mg/L), slightly hard water (150 
mg/L<TH<300 mg/L), hard water (300 mg/L<TH<450 mg/L), and very hard water (TH>450 
mg/L). We compared the TH values with those reported for the major rivers globally and the 
standard values derived from the World Health Organization (WHO, 2011), revealing the 
suitability of the water samples as drinking water. 


3 Results and discussion 


3.1 Hydrochemical characteristics and ion concentration 


The pH, TDS, and ion concentration of the water samples are presented in Table 1. The pH values 
varied from 7.42 to 8.48, with a mean value of 8.00. The TDS ranged from 21.07 to 52.03 mg/L, 
with a mean value of 36.23 mg/L, suggesting an alkaline characteristic and low salinity. The mean 
ion concentration was ranked as follows: Ca?*>SO4? >Na*>NO3 >Mg*'>K*>CI >F >NHs*>Li'*. 
Ca% alone contributed 42% of the total ion concentration, and SOz2”, Na‘, and NO3 accounted for 
41% of the total ion concentration. The cation concentration was ranked as Ca**>Mg?*>Na*>K‘*, 
which indicated that it is the same as the standard crustal material (Ca**>Na‘*>Mg?*>K’*). 


Table 1 Hydrochemical characteristic and ion concentration of river water 


Scien. dene Gi . Cl NO; SO? Lit Nat NH K Mg Ca. 
(mg/L) (mg/L) (mg/L) (mg/L) (mg/L) (mg/L) (mg/L) (mg/L) (mg/L) (mg/L) (mg/L) 

Mean 8.08 40.20 0.06 180 219 532 0.00 3.23 0.08 1.28 1.86 1048 

© SD 013 695 001 042 039 140 0.00 0.84 0.04 0.09 030 174 
Spring Max 831 5203 0.08 253 293 722 000 457 012 143 226 13.53 
Min 7.94 30.51 0.04 126 1.77 291 0.00 192 002 115 137 818 

Mean 7.74 2826 0.05 043 117 283 0.00 221 0.05 O88 118 74l 

SD 0.29 1045 001 027 058 135 0.00 133 0.04 022 042 2.80 
Summer Max 818 49.03 0.07 0.83 229 545 0.00 482 0.13 131 202 12.99 
Min 742 21.07 0.05 0.18 060 1.61 0.00 1.22 0.02 0.73 089 5.66 

Mean 7.94 34.99 081 048 132 430 0.00 269 0.04 1.05 149 9.36 
ummm SD Olt 582 083 013 0.28 075 000 075 002 012 013 1.59 
Max 810 4635 184 072 181 536 001 416 006 127 172 1250 

Min 7.71 3037 0.06 033 0.99 3.09 0.00 201 001 094 135 824 

Mean 822 41.50 0.09 0.95 171 619 0.00 363 016 1.09 1.73 11.06 

l SD 0.19 5.51 0.02 063 0.05 272 0.00 148 0.06 017 010 118 
Winter Max 848 5188 013 217 180 1155 000 654 024 141 183 1329 
Min 8.00 35.63 0.07 040 1.64 3.94 0.00 242 0.09 093 1.54 9.88 


Mean 8.00 36.23 0.25 0.91 160 466 0.00 294 0.09 1.08 157 9.58 
Whole SD 014 7.03 021 035 030 148 0.00 108 003 014 0.22 41.78 
year Max 848 52.03 184 2.53 293 1155 0.01 654 0.24 143 226 13.53 
Min 742 21.07 0.04 0.18 060 161 0.00 1.22 001 0.73 0.89 5.66 


Note: TDS, total dissolved solids; Max, maximum; Min, minimum; SD, standard deviation. 
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Moreover, the anion concentration was ranked as SO4>NO3>CI, which represented a 
difference compared to the standard seawater (CI >SO,?>NO; ). In addition, the order of major 
anions and cations remained the same in different seasons, while the proportion of anions and 
cations varied. Specifically, Ca2*, SOs", Na*, and NO3 accounted for 81% of the total ion 
concentration in spring, the corresponding proportions in summer, autumn, and winter were 84%, 
83%, and 85%, respectively. The result indicated that there are differences in the hydrochemical 
characteristics of river water in different seasons in the study area. 

As depicted in Figure 2, the cations in the water samples were mainly concentrated in the lower 
right corner, indicating a dominance of Ca’* in all seasons. Ca** accounted for 63% of the total 
cations. The proportion of Na* and K* in the total cations reached 26%, while the mean 
proportion of Mg” in the total cations was 10%. In contrast, the anions were characterized by the 
dominance of S04% (63%). NO3 and Cl accounted for 22% and 12% of the total anions, 
respectively. In general, the composition of cations and anions in most of the water samples 
exhibited minute differences across different seasons, demonstrating that the source of these ions 
was natural. 


(a) (b) 0,100 


Spring 
Summer 
Autumn 
Winter 


100 


0 25 50 75 100 
Ca” (%) S0, (%) 


Fig.2 Composition of cations (a) and anions (b) of river water 


3.2 Spatiotemporal variation of the concentration of anions and cations 


As shown in Figure 3, in spring, except for NH4*, the maximum mean concentration of all ions 
occurred at 484 m, while the minimum mean concentration of all ions occurred at 879 m. The 
maximum and minimum mean concentration of NH.* was observed at 537 and 484 m, 
respectively. In summer, except for NH4*, the maximum mean concentration of all ions was 
observed at 879 m, while the minimum mean concentration of all ions mainly occurred at 484 m. 
In autumn, the maximum mean concentration of all ions occurred at 879 m, while the minimum 
mean concentration of all ions mainly occurred at 484 m. In winter, the maximum mean 
concentration of Ca?*, Nat, Mg?*, K+, NHa*, SOu?, Cl, and F- occurred at 879 m, while the 
minimum mean concentration of these ions was mainly observed at 484 m. In general, there were 
significant seasonal and spatial variations in the maximum mean concentration of most ions. 
Furthermore, the minimum mean concentration of most ions was found at 879 m in spring and 
484 m in other seasons. 

The temporal variations of the concentration of anions and cations of river water are illustrated 
in Figure 4. The concentration of CI, NO3, SO4?-, Nat, Mg?*, and Ca’ was characterized by 
significant seasonal variations, with the maximum values observed in winter and the minimum 
values observed in late spring and summer. These ion concentrations varied from 0.35 to 12.51 
mg/L, with a mean value of 3.07 mg/L. Moreover, the ion concentrations were higher in winter 
and early spring than in other seasons. This could be attributed to snowfall accumulation and low 
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runoff in these seasons. The ion concentrations decreased significantly in late spring and summer. 
This result indicated that snowmelt will recharge the river in spring, while precipitation will 
replenish the river in summer. Thus, the major ions in river water were diluted and the ion 
concentration was reduced. 
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Fig. 3 Spatial variations of the concentration of anions (a-d) and cations (e-h) of river water in different seasons 


Concentration of anions (mg/L) 


8.00 


7.00 


6.00 


5.00 


4.00 


3.00 


2.00 


1.00 


q 
= 


15.00 NF 
(b) Cations . 
—_A— NH,* 
—v— K? 

12.00 + A —— Mg” 
T / N¢ —e Ca* d 
=) a 
Š y” 
2 9,00 
8 $ / 
E \ o—¢ 
o Yo 
cas ® 
S / 
2 6.00 + 
3S 
E 
2 

e, 

Q S y 
iS) 3.00 e—e OK eo 4 


Fig. 4 Monthly variations of the concentration of anions (a) and cations (b) of river water 
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3.3 Correlations among pH, electrical conductivity (EC), and the concentration of anions 
and cations 


Correlation analysis is a useful tool for examining the degree of dependence among different 
variables (Anshumali and Ramanathan, 2007). In this study, correlation analysis was performed to 
determine the relationships among pH, EC, and the concentration of anions and cations (Table 2). 
The positive correlation between pH and SO,? indicated that river water is alkaline. The EC 
exhibited a significant positive relationship with all ions except for F~. Moreover, the significant 
correlations of Ca™ with NO3, SO4?, Cl, Na*, K*, and Mg?* suggested that these ions come 
from a common source. No significant correlation was observed between F` and other ions, 
indicating that F` originated from different sources. The presence of Ca? and Mg?’ was primarily 
attributed to the dissolution of carbonates and evaporates. Similarly, the occurrence of SO?” and 
Cl was predominantly linked to the dissolution of evaporates. 


Table 2 Correlation coefficients of pH, electrical conductivity (EC), and the concentration of anions and cations 
of river water 


pH EC F Cl NO; = SO? Li’ Na’ NH,’ K* Mg” Ca** 


pH 1.000 
EC 0.070 1.000 

F 0.059 0.042 1.000 

cl 0.115 0.637" -0.052 1.000 

NO; 0.056 0.792"* -0.049 0.590* 1.000 

SO, 0.203" 0.831 0.104 0.711 0.543" 1.000 

Lit -0.017 0.005 -0.020 -0.048 -0.051 0.013 1.000 

Na‘ 0.145 0.883" 0.041 0.696" 0.593* 0.941" -0.011 1.000 

NH;* 0.055 0.233" -0.031 0.340% 0.013 0.398 -0.061 0.388" 1.000 

K* -0.059 0.761" 0.013 0.657" 0.716"" 0.604 -0.016 0.630% 0.063 1.000 

Mg” 0.093 0.934" 0.061 0.611" 0.866" 0.742" -0.016 0.755"" 0.112 0.800" 1.000 
Ca? 0.199" 0.947" 0.046  0.536™ 0.754"* 0.742* -0.007 0.782" 0.142 0.654" 0.892"" 1.000 


Note: * indicates statistically significant correlation at P<0.01 level; * indicates statistically significant correlation at P<0.05 
level. 


In the PCA, three parameters with eigenvalues greater than 0.600 explained 98% of the total 
variance in the water samples (Table 3). The PC1 (principal component 1; 69%) presented a high 
load of most ions, except for F` and NH4*, indicating a contribution from crustal sources 
(carbonates, evaporates, and sulfates). This was also confirmed by the positive correlations among 
Ca**, Li*, Nat, K*, Mg?*, Cl, SO4*, and NO3- (Table 2). Furthermore, the PC2 (principal 
component 1) accounted for 15% of the total variance, with a high load of NH4*. NH4* was 
strongly correlated with Cl-, SO4?, and Na‘, indicating the same source for these ions. The PC3 
(principal component 3) explained 13% of the total variance, with a high load of F~. The result 
suggested that F~ is mainly related to anthropogenic disturbances. 


3.4 Sources of major ions in river water 


To identify the possible sources of major ions, we conducted the enrichment factor analysis. 
Generally, the EF less than 1.00 indicates dilution compared to the reference material, while the 
EF greater than 1.00 indicates enrichment compared to the reference material (Zhang et al., 2007). 
Furthermore, when the EF is equal to 1.00, it indicates that there is no dilution or enrichment 
relative to the reference material. As depicted in Table 4, the EFsoi value of Na” was less than 
1.00, indicating that the concentration of Na* is primarily influenced by sea. The EFsoi value of 
K* was lower than 1.00, while the EFsea value of K* was greater than 15.00. The result suggested 
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Table 3 Principal component analysis of major ions of river water 


Variable PC1 PC2 PC3 

F —0.077 0.543 0.832 

Cl 0.929 0.104 —0.346 

NO; 0.902 —0.413 —0.018 

S0,” 0.968 0.229 —0.081 

Li? 0.799 0.023 0.484 

Na‘ 0.977 0.067 0.195 

NH,* —0.220 0.771 —0.568 

K* 0.979 0.151 0.002 

Mg” 0.853 0.482 0.191 

Ca?” 0.993 0.072 —0.081 
Explained variance (%) 68.958 15.427 13.334 
Cumulative variance (%) 68.958 84.386 97.719 


Note: PC1, PC2, and PC3 represent principal component 1, principal component 2, and principal component 3, respectively. 


Table 4 Enrichment factor (EF) values of major ions of river water 


EF git EF sea 
Season aS MUMM 
Nat K* Mg” CI NO; SO? K* Mg” Ca? F Cr NO; SO,” 
Spring 0.54 0.24 0.32 55.25 99.15 26.97 18.21 2.55 74.04 28.34 0.48 456.94 13.62 
Summer 0.52 0.24 0.28 18.56 74.78 20.29 18.23 2.35 76.33 34.77 0.17 355.33 10.56 
Autumn 0.50 0.22 0.28 16.67 67.21 24.46 17.94 2.44 79.31 435.29 0.16 331.84 13.23 
Winter 0.58 0.20 0.28 27.72 29.77 29.77 13.80 2.09 69.36 34.55 0.23 317.15 14.08 


Whole year 0.54 0.22 0.29 30.78 79.22 25.88 16.79 2.34 74.23 124.49 0.27 366.04 13.10 


Note: EFsoi represents the enrichment factor for soil, and EFsea represents the enrichment factor for sea. 


that the concentration of K* is regulated by both sea and soil. The EFsoi value of Mg?* was lower 
than 1.00, while its EFsea value was higher than 1.00, suggesting that Mg”* is derived from both 
soil and sea. The EFsoi value of Ca** ranged from 69.36 to 79.31, indicating that Ca?* 
predominantly comes from soil. In addition, F- was affected by soil, as evidenced by its EFsoi1 
value. The EF sea value of Cl was below 1.00, while its EF oi value was higher than 10.00. This 
result demonstrated that the concentration of Cl is regulated by both sea and soil. Furthermore, 
NO3° was mainly affected by soil, with both EF,oi value and EFsea value of NO3 above 1.00. 
SO, was mainly affected by soil, with an EFsoi exceeding 1.00. In addition, for all ions, 
differences were observed among different seasons. 


3.5 Factors influencing hydrochemical characteristics of rivers 


The hydrochemical characteristics of rivers can be influenced by various anthropogenic activities, 
such as increased grazing intensity, which has been known to modify water chemistry in recent 
years (Goenster-Jordan et al., 2021). Additionally, it has been reported that the growth of tourism 
and local agricultural practices have the potential to alter natural water chemistry (Negrel et al., 
1993). Therefore, in order to assess the impact of rock weathering and human activities on water 
quality, we calculated the ratios of NOx, SO4%, and CT to Na* according to Gaillardet et al. 
(1999a, b). Figure 5 illustrates that no significant correlation was observed between the ratio of 
NO3 to Na* (NO3/Na’*) and the ratio of S04% to Na* (S04%/Na*) in different seasons, suggesting 
that the sources of these ions were different. Most of the NO3/Na* values were below 0.8, 
indicating a minimal contribution from human activities. 

Furthermore, we drew the Gibbs diagram based on the ratio of Na* to Na* and Ca?* (Na*/(Na*+ 
Ca**)). The TDS content was employed to identify the predominant factors influencing 
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hydrochemical characteristics of rivers such as precipitation dominance, weathering dominance, 
or evaporation-crystallization dominance. As shown in Figure 6, most of the water samples had 
moderate TDS and low Na*/(Na*+Ca’*) values, indicating that rock weathering was the leading 
process regulated hydrochemical characteristics of rivers. To distinguish the main rocks (i.e., 
carbonates, evaporates, and sulfates), we calculated the ratio of Ca** to Na* (Ca**/Na*) and the 
ratio of Mg% to Na” (Mg?*/Na*). The majority of the water samples were found to lie between 
carbonate dissolution and silicate weathering, indicating that both processes contribute 


significantly to the hydrochemical characteristics of rivers. 
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Fig. 5 (a), relationship between the ratio of NO37 to Na* (NO3"/Na‘*) and the ratio of S04% to Na* (S04%/Na*); 
(b), relationship between the ratio of CI- to Na* (CI-/Na*) and NO3-~/Na* 
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Fig. 6 Influencing factors for hydrochemical characteristics showed by Gibbs diagram. (a), Gibbs diagram of 
the ratio of Na* to Na* and Ca% (Nat/(Na™+Ca*)) to total dissolved solids (TDS); (b), relationship between the 
ratio of Ca? to Na* (Ca*/Na*) and the ratio of Mg** to Na* (Mg?*/Na‘). 


The main ions of water river in the western region of the Altay Mountains were Ca’, Mg”, 
and SO, (Fig. 2). However, significant differences in ion concentration were observed as the 
elevation gradient changed (Fig. 3). At low elevations (<500 m), Ca** and Cl were the dominant 
ions; while, at high elevations (from 500 to 900 m), Ca? and SO4% were the dominant ions. A 
related study observed that the vertical distribution of major ions can provide insights into the 
vertical distribution of rock formations (Li et al., 2018). The maximum values of K*, NH4*, F-, 
Li‘, and Mg** remained consistent across different seasons, suggesting that the concentrations of 
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these ions were minimally affected by precipitation and meltwater recharge. The maximum values 
of Cl’, NO37, S047, Ca**, and Na’ first decreased and then increased from spring to winter. The 
result indicated that the concentrations of Cl, NO7, SO4%, Ca**, and Nat are affected by the 
recharge effects of precipitation and meltwater. However, the recharge effects of precipitation and 
meltwater on K*, NH4*, F-, Li*, and Mg?* were minimal. 


3.6 Comparison of hydrochemical characteristics between the study area and other mountains 


The hydrochemical characteristics of the study area are consistent with those of other mountains 
(Table 5). For instance, the pH values of the Ertix River and Ulungur River were also alkaline, 
with a small range of variation. These rivers contained abundant Ca?*, Na*, HCO3,, and SO,” , 
which are the most common cations and anions (Liu et al., 2021). Additionally, the average TDS 
in the Ertix River and Ulungu River is reported to be significantly higher than that in the western 
region of the Altay Mountains, which may be attributed to the comprehensive impact of 
geological conditions, supply sources, climatic factors, and human activities (Zabaleta et al., 
2007). The hydrochemical characteristics of the Burqin River and Haba River were related to the 
weathering of carbonates and silicates. In contrast, it has been reported that the ion composition 
of the Ertix River and Ulungu River is influenced by various factors, including carbonates 
dissolution, silicates weathering, and evaporates dissolution (Liu et al., 2021). 


Table 5 Comparison of hydrochemical characteristics between the western region of the Altay Mountains and 
other mountains 


TDS Clr NO; SO, HCO; Na? K* Mg* Ca” 


inch ial PH (mg/L)(mg/L)(mg/L)(mg/L)(mg/L)(mg/L)(mg/L)(mg/L) (mg/L) _ Reference 
Bisa Riverand 8.00 36.20 0.91 160 4.66 - 2.94 1.08 1.57 9.58 This study 
Altay Haba River 
Mountains ie Ri 
Bei Rivet and 7.80 126.00 7.06 - 35.00 65.70 12.60 1.50 4.26 32.40 Liu etal. (2021) 
Ulungur River 


Source region of SYF $ 20 291.40 3.00 3.50 45.00200.00 5.00 1.44 11.00 35.00 Ma etal. (2019) 


f Darya River 

Tianshan 

Mountains Kax River 8.00 239.20 1.70 - 22.55120.80 4.10 1.67 5.83 39.38 Feng et al. (2022) 
Bortala River 8.50 246.00 14.30 - 67.50157.00 21.40 2.32 7.83 52.80 Liu ctal. (2021) 

Qilian Hulugou River - = 16.60 2.67 360.05 - 82.67 2.29 33.70 88.21 Li etal. (2014) 

Mountains Shule River 8.50 160.10 8.20 - 23.70 88.70 12.60 0.90 10.40 14.80 Qu etal. (2019) 

Central Source region of 7.60 778.00233.70 1.30 114.90 188.50 157.70 5.50 22.90 53.40 Jiang et al. (2015) 


Tibetan Yangtze River 
Plateau Zhajia Tsangpo River 9.10 947.00 394.81 2.62 164.11 55.20245.00 12.20 27.30 45.50 Qu etal. (2019) 


Southern and Yellow River 8.40 609.50 22.90 - 49.40 389.50 33.30 2.10 33.70 75.80 Qu etal. (2019) 
T Heihe River 8.70 462.00 8.20 2.20 80.40263.30 15.70 1.70 30.40 60.00 Qu etal. (2019) 
parts of the 

Tibetan Nujiang River 8.40 141.00 5.00 0.00 31.00 66.00 3.00 1.00 7.00 24.00 Huang etal. (2009) 


Plateau Yarlung Zangbo River 8.80 117.30 2.40 1.00 37.40 74.10 5.70 1.00 5.10 27.90 Quetal. (2017) 


Note: -, no data. 


The ion concentration of river water in the western region of the Altay Mountains was far lower 
than that in other mountains, which is attributed to the different and unique lithology and 
geographic conditions of these mountains (Table 5). For instance, the mean concentration of Cl- 
in the study area was 0.91 mg/L, which is lower than that reported in other mountains. The TDS 
of the study area was also relatively low, due to lower Ca?* and SO4%. The mean concentrations 
of Ca?* and SO,?- were 9.58 and 4.66 mg/L, respectively, which are more than 17 times lower 
than those in other rivers. In addition, the high concentrations of Ca?*, HCO3°, and SO4% in the 
rivers of the Tianshan Mountains demonstrated that hydrochemical characteristics are mainly 
affected by carbonate dissolution and silicate weathering (Ma et al., 2019; Liu et al., 2021; Feng 
and Yang, 2022). High concentrations of Ca?*, Mg**, and HCO; in rivers in southern and 
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northern parts of the Tibetan Plateau were also affected by carbonate dissolution and silicate 
weathering (Huang et al., 2009; Qu et al., 2017, 2019). However, high concentrations of Na*, Cl, 
and SO4% in rivers of central Tibetan Plateau were controlled by evaporates dissolution (Galy and 
France-Lanord, 1999). Furthermore, high concentrations of Ca?*, Mg?*, and SO, in the Hulugou 
River of the Qilian Mountains were affected by the weathering and dissolution of dolomite (Li et 
al., 2014). However, it is reported that in the Shule River of the Qilian Mountains, the ion 
composition (Ca?*, Mgt, and HCO;)) is mainly influenced by the weathering of carbonate and 
silicate (Qu et al., 2019). 


3.7 Impacts on the suitability of river water for drinking and irrigation 


Glaciers and snowmelt in the Altay Mountains are important water resources for residents and 
economic development of this region (Wang et al., 2015). However, due to the high sensitivity of 
small glaciers in the Altay Mountains, the lack of water may threaten this region. Therefore, 
considering the impact of water scarcity on human and ecosystem health, we conducted an 
assessment of river water quality. In accordance with the guidelines of the World Health 
Organization, we found that the TDS and ion concentrations in the study area were considerably 
below the maximum limits for drinking water (WHO, 2011). The water quality of this study area 
is similar to that of the Indus River and Ganges River—Yarlung Zangbo River (Qu et al., 2019). 
Furthermore, all water samples were very soft water, indicating that the water is safe to drink 
(Fig. 7). 


10000 
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100 S10 


S11 
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10 100 1000 10000 
TH (mg/L) 
Fig. 7 River water quality at different sampling sites in the western region of the Altay Mountains. TH represent 
the total hardness. S1 to S12 represent 12 sampling sites, and the locations of these sites are depicted in Figure 1. 


We also assessed the suitability of river water for irrigation (Table 6). The mean SSP value in 
the region was 18% (+4%) (Table 6), indicating excellent irrigation quality. Moreover, compared 
to other sampling sites, the SSP of S1 (Site 1) was relatively high, which may be related to human 
activities in the downstream region. In this study, most of the SPP values (81%) achieved 
excellent irrigation water quality. Similarly, the mean SAR value was 1.2 (+0.5), and the irrigation 
water quality at all sampling sites was excellent (Table 6). Overall, river water in the western 
region of the Altay Mountains is safe for drinking and irrigation. 

We evaluated the suitability of river water for drinking and irrigation based on the TDS, TH, 
SSP, and SAR of the water samples. Future studies could also focus on comprehensive long-term 
field observations of water chemistry and isotopes in glacial areas. This will facilitate the 
systematic assessment of water quality and hydrochemical characteristics in headwater 
catchments, particularly in data-poor regions, such as the Altay Mountains. 
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Table 6 Soluble sodium percentage (SSP) and sodium adsorption ratio (SAR) of river water in the western 
region of the Altay Mountains 


Sampling site SSP (%) SAR Sampling site SSP (%) SAR 


Sl 21+5 7 S7 18+6 1.1 
S2 18+2 1.1 S8 1743 1.0 
S3 17+2 1.1 S9 203 1.3 
S4 19+3 1.2 S10 163 0.9 
S5 16+4 0.9 S11 18+6 1.0 
S6 2247 1.7 S12 183 1.2 


Note: Mean+SD. 


4 Conclusions 


In this study, the mean ion concentration was ranked as follows: Ca?*>SO4* >Na*>NO3>Mg?*> 
K*>Cl->F >NH4‘>Li*. Among these ions, Ca?*, SO4?-, Na*, and NO37 accounted for 83% of the 
total ion concentration. The spatial variation of ion concentration exhibited greater stability in 
summer, autumn, and winter. Results of Gibbs diagram showed that river water chemistry was 
primarily influenced by carbonate dissolution and silicate weathering. There were differences in 
hydrochemical characteristics at different elevations, with low elevations dominated by Ca?* and 
Cl, while high elevations dominated by Ca** and SO4%. In addition, the recharge effects of 
precipitation and snowmelt primarily affected the concentrations of CI, NOs, SOq?, Ca?*, and 
Na*. The suitability evaluation of river water revealed that river water is safe for drinking and 
irrigation. Nonetheless, this study is only a preliminary investigation, and follow-up researches 
may focus on further strengthening the field monitoring of water chemistry and isotopes in the 
glacial region of the Altay Mountains. 
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